Frictional Melting in Hydrothermal Fluid-Rich Faults: Field and Experimental Evidence From the Bolfín Fault Zone (Chile).

Tectonic pseudotachylytes are thought to be unique to certain water-deficient seismogenic environments and their presence is considered to be rare in the geological record. Here, we present field and experimental evidence that frictional melting can occur in hydrothermal fluid-rich faults hosted in the continental crust. Pseudotachylytes were found in the >40 km-long Bolfín Fault Zone of the Atacama Fault System, within two ca. 1 m-thick (ultra)cataclastic strands hosted in a damage-zone made of chlorite-epidote-rich hydrothermally altered tonalite. This alteration state indicates that hydrothermal fluids were active during the fault development. Pseudotachylytes, characterized by presenting amygdales, cut and are cut by chlorite-, epidote- and calcite-bearing veins. In turn, crosscutting relationship with the hydrothermal veins indicates pseudotachylytes were formed during this period of fluid activity. Rotary shear experiments conducted on bare surfaces of hydrothermally altered rocks at seismic slip velocities (3 m s-1) resulted in the production of vesiculated pseudotachylytes both at dry and water-pressurized conditions, with melt lubrication as the primary mechanism for fault dynamic weakening. The presented evidence challenges the common hypothesis that pseudotachylytes are limited to fluid-deficient environments, and gives insights into the ancient seismic activity of the system. Both field observations and experimental evidence, indicate that pseudotachylytes may easily be produced in hydrothermal environments, and could be a common co-seismic fault product. Consequently, melt lubrication could be considered one of the most efficient seismic dynamic weakening mechanisms in crystalline basement rocks of the continental crust.

Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone.

In subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible for the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the Papaku thrust (Hikurangi subduction zone, New Zealand). Here, we show that at seismic velocity, shear-induced dilatancy is followed by pressurisation of fluids. The thermal and mechanical pressurisation of fluids, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.

Mechanical behaviour of fluid-lubricated faults.

Fluids are pervasive in fault zones cutting the Earth's crust; however, the effect of fluid viscosity on fault mechanics is mainly conjectured by theoretical models. We present friction experiments performed on both dry and fluid-permeated silicate and carbonate bearing-rocks, at normal effective stresses up to 20 MPa, with a slip-rate ranging between 10 µm/s and 1 m/s. Four different fluid viscosities were tested. We show that both static and dynamic friction coefficients decrease with viscosity and that dynamic friction depends on the dimensionless Sommerfeld number (S) as predicted by the elastohydrodynamic-lubrication theory (EHD).Under favourable conditions (depending on the fluid viscosity (η), co-seismic slip-rate (V), fault geometry (L/H02) and earthquake nucleation depth (∝σeff)), EHD might be an effective weakening mechanism during natural and induced earthquakes. However, at seismic slip-rate, the slip weakening distance (Dc) increases markedly for a range of fluid viscosities expected in the Earth, potentially favouring slow-slip rather than rupture propagation for small to moderate earthquakes.

Subseismic to Seismic Slip in Smectite Clay Nanofoliation.

Smectite clays are the main constituent of slipping zones found in subduction zone faults at shallow depth (e.g., <1-km depth in the Japan Trench) and in the decollements of large landslides (e.g., 1963 landslide, Vajont, Italy). Therefore, deformation processes in smectite clays may control the mechanical behavior from slow creep to fast accelerations and slip during earthquakes and landslides. Here, we use (1) laboratory experiments to investigate the mechanical behavior of partly water-saturated smectite-rich gouges sheared from subseismic to seismic slip rates V and (2) nanoscale microscopy to study the gouge fabric. At all slip rates, deformation localizes in volumes of the gouge layer that contain a "nanofoliation" consisting of anastomosing smectite crystals. "Seismic" nanofoliations produced at V = 0.01, 0.1, and 1.3 m/s are similar to "subseismic" nanofoliations obtained at V = 10-5 m/s. This similarity suggests that frictional slip along water-lubricated smectite grain boundaries and basal planes may occur from subseismic to seismic slip rates in natural smectite-rich faults. Thus, if water is available along smectite grain boundaries and basal planes, nanofoliations can develop from slow to fast slip rates. Still, when nanofoliations are found highly localized in a volume, they can be diagnostic of slip that occurred at rates equal or larger than 0.01 m/s. In such a case, they could be markers of past seismic events when found in natural fault rocks.

An empirically based steady state friction law and implications for fault stability.

Empirically based rate-and-state friction laws (RSFLs) have been proposed to model the dependence of friction forces with slip and time. The relevance of the RSFL for earthquake mechanics is that few constitutive parameters define critical conditions for fault stability (i.e., critical stiffness and frictional fault behavior). However, the RSFLs were determined from experiments conducted at subseismic slip rates (V < 1 cm/s), and their extrapolation to earthquake deformation conditions (V > 0.1 m/s) remains questionable on the basis of the experimental evidence of (1) large dynamic weakening and (2) activation of particular fault lubrication processes at seismic slip rates. Here we propose a modified RSFL (MFL) based on the review of a large published and unpublished data set of rock friction experiments performed with different testing machines. The MFL, valid at steady state conditions from subseismic to seismic slip rates (0.1 µm/s < V < 3 m/s), describes the initiation of a substantial velocity weakening in the 1-20 cm/s range resulting in a critical stiffness increase that creates a peak of potential instability in that velocity regime. The MFL leads to a new definition of fault frictional stability with implications for slip event styles and relevance for models of seismic rupture nucleation, propagation, and arrest.

G: Fracture energy, friction and dissipation in earthquakes.

Recent estimates of fracture energy G' in earthquakes show a power-law dependence with slip u which can be summarized as G' ∝ ua where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to Gf : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G' and Gf should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s2). The resulting fracture energy Gf is similar to the seismological estimates, with Gf and G' being comparable over most of the slip range. However, Gf appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G' appears to increase further and surpasses Gf at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.

Dynamic weakening of serpentinite gouges and bare surfaces at seismic slip rates.

To investigate differences in the frictional behavior between initially bare rock surfaces of serpentinite and powdered serpentinite ("gouge") at subseismic to seismic slip rates, we conducted single-velocity step and multiple-velocity step friction experiments on an antigorite-rich and lizardite-rich serpentinite at slip rates (V) from 0.003 m/s to 6.5 m/s, sliding displacements up to 1.6 m, and normal stresses (σn ) up to 22 MPa for gouge and 97 MPa for bare surfaces. Nominal steady state friction values (µnss) in gouge at V = 1 m/s are larger than in bare surfaces for all σn tested and demonstrate a strong σn dependence; µnss decreased from 0.51 at 4.0 MPa to 0.39 at 22.4 MPa. Conversely, µnss values for bare surfaces remained ∼0.1 with increasing σn and V. Additionally, the velocity at the onset of frictional weakening and the amount of slip prior to weakening were orders of magnitude larger in gouge than in bare surfaces. Extrapolation of the normal stress dependence for µnss suggests that the behavior of antigorite gouge approaches that of bare surfaces at σn ≥ 60 MPa. X-ray diffraction revealed dehydration reaction products in samples that frictionally weakened. Microstructural analysis revealed highly localized slip zones with melt-like textures in some cases gouge experiments and in all bare surfaces experiments for V ≥ 1 m/s. One-dimensional thermal modeling indicates that flash heating causes frictional weakening in both bare surfaces and gouge. Friction values for gouge decrease at higher velocities and after longer displacements than bare surfaces because strain is more distributed. KEY POINTS: Gouge friction approaches that of bare surfaces at high normal stressDehydration reactions and bulk melting in serpentinite in < 1 m of slipFlash heating causes dynamic frictional weakening in gouge and bare surfaces.

Fault lubrication during earthquakes.

The determination of rock friction at seismic slip rates (about 1 m s(-1)) is of paramount importance in earthquake mechanics, as fault friction controls the stress drop, the mechanical work and the frictional heat generated during slip. Given the difficulty in determining friction by seismological methods, elucidating constraints are derived from experimental studies. Here we review a large set of published and unpublished experiments (∼300) performed in rotary shear apparatus at slip rates of 0.1-2.6 m s(-1). The experiments indicate a significant decrease in friction (of up to one order of magnitude), which we term fault lubrication, both for cohesive (silicate-built, quartz-built and carbonate-built) rocks and non-cohesive rocks (clay-rich, anhydrite, gypsum and dolomite gouges) typical of crustal seismogenic sources. The available mechanical work and the associated temperature rise in the slipping zone trigger a number of physicochemical processes (gelification, decarbonation and dehydration reactions, melting and so on) whose products are responsible for fault lubrication. The similarity between (1) experimental and natural fault products and (2) mechanical work measures resulting from these laboratory experiments and seismological estimates suggests that it is reasonable to extrapolate experimental data to conditions typical of earthquake nucleation depths (7-15 km). It seems that faults are lubricated during earthquakes, irrespective of the fault rock composition and of the specific weakening mechanism involved.